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We study the charged non-relativistic Bose gas interacting 
with a constant magnetic field but which is otherwise free. 
The notion of Bose-Einstein condensation for the three di- 
mensional case is clarified, and we show that although there 
is no condensation in the sense of a phase transition, there is 
still a maximum in the specific heat which can be used to de- 
fine a critical temperature. Although the absence of a phase 
transition persists for all values of the magnetic field, we show 
how as the magnetic field is reduced the curves for the spe- 
cific heat approach the free field curve. For large values of 
the magnetic field we show that the gas undergoes a "dimen- 
sional reduction" and behaves effectively as a one-dimensional 
gas except at very high temperatures. These general features 
persist for other spatial dimensions D and we show results 
for D — 5. Finally we examine the magnetization and the 
Meissner-Ochsenfeld effect. 



The thermodynamical properties of a charged non- 
relativistic Bose gas in a constant magnetic field were 
correctly evaluated by Schafroth Q . He showed that the 
presence of a constant magnetic field destroyed the famil- 
iar Bose-Einstein condensation exhibited by the free gas. 
Furthermore, Schafroth demonstrated that the magne- 
tized gas exhibited the Meissner-Ochsenfeld effect which 
set in at a scale related to the condensation temperature 
of the free Bose gas. It is important to emphasize what 
is meant by the absence of Bose-Einstein condensation 
for the magnetized Bose gas. Unlike the situation for the 
free Bose gas, there is no phase transition (in the infinite 
volume limit). However this does not imply that there 
is no build-up of particles in the ground state. (This 
point has been emphasized in Ref. |2|3j.) Subsequent to 
the work of Schafroth, May |Q looked at the magnetized 
Bose gas in a general spatial dimension D and showed 
that for D > 5 there was a phase transition like that in 
the free Bose gas. Thus the study of spatial dimensions 
other than 3 can be used to illustrate general features of 
phase transitions in magnetized systems. 

We will begin with the expression for the total ther- 
modynamic potential as found using the effective action 
method. (For a review see @j.) This is 
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The first term in Vt is the "classical" part which accounts 
for a possible condensate described by and the second 
term is the standard expression for the thermodynamic 
potential from statistical mechanics. Here (3 = (fcT) -1 , /i 
is the chemical potential, and E n are the energy levels of 
the system. Q, in ([!]) can be evaluated from a knowledge 
of the energy levels for a charged particle in a constant 
magnetic field || . Although in D spatial dimensions the 
magnetic field may in general have more than one inde- 
pendent component, we will only consider the presence 
of a single component here. 

Whether or not there is a non-zero condensate ^ ^ 
is determined by the solution to 

_ J-D 2 ^ - e/xf = . (2) 
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This is obtained by extremizing (Q) with respect to , F. It 
is easy to show that = if /i < fx c where e/i c = E gives 
a critical value for the chemical potential determined by 
the lowest energy eigenvalue. For the case of a constant 
magnetic field of strength B we have Eq = lu/2 where u = 
eB/m. In this case there is no non-trivial condensate. 
(See Ref. || for a general treatment.) If it is possible for 
fi to reach the critical value fi c , then we have \I/ oc /o(x) 
where /o(x) is the eigenfunction corresponding to the 
lowest energy. 

Whether or not \x can reach yit c is determined by the 
total charge Q = — |p. If we set the condensate # = 0, 
so that the first term in f2 makes no contribution to the 
charge, and if we can solve for /j, in terms of Q, then we 
must have /i < fi c . However if it is not possible to do this 
then "J is necessarily non-zero and we have a condensate. 
The nature of the solutions for [i changes with temper- 
ature, and it is this behaviour which can give rise to a 
critical temperature signalling a phase transition. The 
case of a condensate ( 1 F ^ 0) and no condensate = 0) 
must be distinguished. From previous work [Q,^ we know 
that when a magnetic field is present the condition for 
# ^ is D > 5. This contrasts with the condition for 
the free gas (B = 0) which requires the spatial dimension 
D > 3 for a non-zero condensate. 

In the physically interesting case of D — 3 the presence 
of a magnetic field destroys Bose-Einstein condensation 
if it is interpreted as a phase transition with a non-zero 
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condensate Q. The presence of a magnetic field, no mat- 
ter how small, results in the charged Bose gas behaving in 
a different manner from the same system when no mag- 
netic field is present. However on physical grounds one 
might expect that for a very small magnetic field the sys- 
tem should behave in almost the same way as it would if 
B = 0. This is one of the issues addressed in the present 
paper. In addition the approach of Ref. || is clarified 
with a particular example. 

For the free Bose gas a key feature of Bose-Einstein 
condensation is the behaviour of the specific heat at the 
critical temperature To. For D = 3 the specific heat has 
a maximum at T = To, and is continuous there with a 
discontinuous first derivative ||. For D = 4 the specific 
heat and its first derivative are continuous at T = To, and 
for D > 5 the specific heat is discontinuous at T = Tq ||. 
In all cases the non-smooth behaviour of the specific heat 
corresponds to a maximum. As already stated the pres- 
ence of a magnetic field suppresses the phase transition if 
D = 3,4. For D > 5 the specific heat for the magnetized 
gas would be expected to show a non-smooth behaviour 
at a critical temperature T c corresponding to the exis- 
tence of a phase transition; however for D = 3, 4 the 
specific heat should be perfectly smooth since no phase 
transition occurs. We will concentrate initially on D = 3. 

By first calculating the internal energy defined by U — 
dfl/(d(3) with n/3,V,B held fixed, and then the specific 
heat at constant volume by Cy = dU / (dT) with V, B 
and the total charge Q held fixed, it is possible to show 
that 
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is the total number of particles present. Here ip(a,b,c) 
denotes a class of sums defined by 
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(We take D = 3 in eqs.(|) and §).) 

It is possible to solve (||) for e where fi = w(l/2 — e) 
and then to calculate Cv ■ The result of this is shown 
m fig. [J For comparison we have also indicated the re- 
sult for the free Bose gas. As expected the specific heat 
is perfectly continuous for B ^ 0. As B is decreased 
(following the curves with smaller values of Xo) although 
Cy remains continuous the peak starts to increase and 



the curve starts to resemble that for the free Bose gas. 
The maximum in the specific heat becomes sharper as B 
is decreased with the temperature of the maximum ap- 
proaching the temperature To for the free Bose gas. Sim- 
ilar results hold if the number of particles in the ground 
state is computed. Thus although for B ^ the sys- 
tem does not have a phase transition characterized by a 
critical temperature and a non-zero condensate "J and is 
therefore very different from the same system with B = 0, 
quantities of direct physical interest, such as the specific 
heat or ground state particle number, are closely approx- 
imated by the free Bose gas result for small enough values 
of B. We will obtain approximate analytical expressions 
for small B later. 
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FIG. 1. This shows the specific heat for the magnetized 
Bose gas in units of kN. The curve labelled free is the result 
for the uncharged Bose gas. The x coordinate in the figure is 
x = huj/(kT) in units of xq — hu/(kTo) with To the critical 
temperature for the free Bose gas. The numbers labelling the 
curves are the values for xo. The smaller numbers indicate 
lower values for the magnetic field. 

It is also of interest to observe what happens as the 
value of the magnetic field is increased. It is clear from 
the curves in fig. |l| with larger values of xo that with 
increasing B the peak in Cy shifts to smaller values of 
x (which corresponds to higher temperatures) and de- 
creases in magnitude. (Although the resolution of the 
graphs do not show it, there is always a maximum in the 
specific heat.) As x approaches 0, corresponding to in- 
finitely high temperatures, the specific heat approaches 
the classical Maxwell-Boltzmann result of 1.5 regardless 
of the magnitude of J5, as would be expected. For large B 
the specific heat actually resembles that for a free Bose 
gas in one spatial dimension. This can be understood 
heuristically by the fact that classically the orbit of a 
charged particle in a constant magnetic field is a spiral 
around the direction specified by the magnetic field. The 
size of the orbit is inversely proportional to the magnetic 
field; thus as B is increased the motion of the particle be- 
comes more and more one-dimensional. This substanti- 
ates the general point of view used in Ref. || in which the 
leading behaviour of thermodynamic quantities is well 
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approximated by studying the lowest energy solution. In 
this case as B is increased the gap between the ground 
and first excited state becomes increasingly large and the 
ground state makes the leading contribution. 

Although for B ^ with D = 3 there is no phase 
transition, Schafroth [Q showed that the magnetized gas 
exhibited the Meissner-Ochsenfeld effect. He found 
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where p is the charge density and xo — to/ (kTo). The ap- 
proximation assumes that T < Tq. Using M = d£l/ (dB), 
the exact expression may be computed to be 



M = -- 
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This expression holds for any temperature. 
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FIG. 2. This shows the ratio of the magnetization com- 
puted using the Schafroth approximation in eq. ^ to the 
exact result. The labels 100,10,1,0.1 on the curves indicate 
the value of xq, the smaller numbers corresponding to lower 
values for the magnetic field. 

The ratio of the Schafroth approximation to the true 
result is plotted in fig. ^. As B is lowered the Schafroth 
approximation becomes better over a wider range of tem- 
perature, and as T — » Tq. As noted by Schafroth M, the 
approximation will break down when T becomes too close 
to To, and this is borne out by the detailed calculation. 
For T ~ To/10 the agreement is within 2%. 

It is possible to derive the Schafroth approximation in 
(0) from the exact result (fy by approximating for the 
class of sums (|E]). A powerful technique for doing this 
makes use of the Mellin-Barnes integral representation 
for the exponential function Jl0|,[| . In Ref. |l(J it was 
shown that for e ~ 1/2 (corresponding to /i ~ 0) there 
was a superdiamagnetic regime where the magnetization 
M oc —B 1 / 2 . (A previous paper which studied the mag- 
netization is Ref. pH , although the superdiamagnetism 
was not found.) The temperature range at which e ~ 1/2 
is close to but slightly above Tq. We have obtained re- 
sults in a temperature range which includes T = Tq. By 
solving fl) for e we find 



e ~ a + 
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which is good if 1 — (x/xo) 1 ^ 2 << x^ 2 . Here a is just a 
constant defined by £#(1/2, a) = with the Hurwitz 
C-function. (Numerically we find a ~ 0.302721829.) As 
xo is decreased the range over which the expansion holds 
becomes increasingly small, but it is always good at x = 
xq (corresponding to T = To). The agreement between 
the approximation (j^) and a numerical evaluation of e is 
very good for small xq. 
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FIG. 3. This shows the specific heat for small values of the 
magnetic field. The labels on the curve indicate the results 
for xo = 0.1, 0.5, 1. For comparison the result for no magnetic 
field is shown in the discontinuous curve. 

Although it might be thought that there was very little 
difference between our result of e ~ 0.3 and the result 
e ~ 0.5 of Ref. |{Tof , in fact the magnetization and other 
thermodynamic quantities are extremely sensitive to the 
value of e. Using our expansion (ra) we find 

\f - . 9 f 67rV a C H (-l/2,o) i/ 2 
2ml C(3/2) 

oC(l/2U 
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for small xq. This shows the —B 1 ^ 2 behaviour found in 
Ref. [0, but with a slightly different numerical factor 
due to the lower temperature range we are using. 

In addition to the magnetization we have obtained an 
approximate expression for the specific heat using (||) 
which is 



C v _ 15 C(5/2) _ 
kN ~ 4 C(3/2) _ 27r 1 /2^(3/2,a)■ 
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at T = T). The first term in (10) which survives in the 
B — > limit may be observed ]8f to be the specific heat 
for the free Bose gas at the critical temperature Tq. The 
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presence of a non-zero magnetic field is seen to lower the 
specific heat from the free field result which is consistent 
with the results shown in fig. 1. It is possible to system- 
atically determine higher order terms in the expansions 
(|s|-|lO|) and increase the accuracy of the approximation. 
This and other details will be presented elsewhere |T^| . 

As an example of a magnetic system where a phase 
transition does occur we will consider the case D = 5. 
When B = the specific heat is discontinuous at the 
critical temperature To. For B ^ the specific heat is 
given by an expression similar to (|^), and is perfectly 
continuous. The results are shown in fig. ^. 

The critical temperature T c for the phase transition 
cannot be computed analytically if B ^ 0; however for 
small values of the magnetic field it is possible to show 
that § 



kT c ~ kT 



COT eB 
5C(5/2) to 
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The critical temperature is lowered from that for the free 
Bose gas when there is a non-zero magnetic field. The ap- 
proximation in (|ll]) is in good agreement with the exact 
result for small B but breaks down as B becomes large 
|l2|] . The results in fig. || show that as B is decreased 
the curve for the specific heat starts to resemble the dis- 
continuous one for the free gas. The critical temperature 
approaches the free field result, which is clear from the 
small B approximation in (|Tl|). The height of the spe- 
cific heat maximum is seen to increase in magnitude and 
the slope of the curve for T > T c becomes steeper as 
B is reduced. On the other hand as the magnetic field 
is increased in magnitude the specific heat undergoes a 
marked change. The results for larger values of the mag- 
netic field are shown in fig. 0. 
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FIG. 4. This shows the specific heat for small values of the 
magnetic field. The labels on the curve indicate the results 
for x = 1,2,3,5,10. 



In conclusion, we have presented the results of a study 
of the ideal charged Bose gas in a constant magnetic field. 
For certain values of the temperature and magnetic field 
we have obtained analytical results and compared them 
with numerical ones. By examining the specific heat it 
was shown how in a strong magnetic field the system 
exhibits a "dimensional reduction" where the system be- 
haves as if the spatial dimension was lowered by 2. We 
have also examined the way in which the free Bose gas 
results are approached for small B. Although there is 
a significant difference between the case B = (which 
has a phase transition) and the B ^ gas (which has 
no phase transition), for small values of B there is very 
little difference between the two specific heat curves. Fi- 
nally for spatial dimensions D > 5 a proper account of 
the condensate is essential for obtaining correct results, 
and this can be done in an efficient manner by using the 
effective action method. 
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For very large values of B the curves for the specific 
heat resemble the result found for the free Bose gas in 
three spatial dimensions. This provides a further and 
striking example of the "dimensional reduction" in the 
presence of a magnetic field. 
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